Andreev and Single Particle Tunneling Spectroscopies in Underdoped Cuprates 
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We study tunneling spectroscopy between a normal metal and underdoped cuprate superconductor 
modeled by a phenomenological theory in which the pseudogap is a precursor to the undoped Mott 
insulator. In the transparent tunneling limit, the spectra show a small energy gap associated with 
Andreev reflection. In the Giaever limit, the spectra show a large energy gap associated with 
single particle tunneling. Our theory semi-quantitatively describes the two gap behavior observed 
in tunneling experiments. 
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It has long been argued that the highly anomalous 
pseudogap phase of underdoped (UD) cuprates holds the 
key to the physics of high T^ superconductors [IHl]- A 
variety of models have been proposed to describe this 
phase. At present a consensus is lacking and the mer- 
its of the different models are being vigorously debated. 
Some models propose the partial truncation of the Fermi 
surface in the pseudogap phase is due to the presence of 
order, which breaks translational symmetry. However in 
the absence of experimental evidence for broken transla- 
tional symmetry [S], models without this property have 
gained traction. These in turn can be divided into two 
classes. One emphasizes the reduction of the supercon- 
ducting (SC) Tc by strong phase fluctuations due to the 
reduced superfluid density in UD cuprates. This allows 
the larger SC energy gaps, which in the d-wave form are 
maximal in the antinodal directions on the Fermi sur- 
face, to remain finite at temperatures T > Tc due to 
preformed Cooper pairs but without off-diagonal long 
order 0. Nernst effect and diamagnetism experiments 
confirm the presence of SC fiuctuations in an extended 
temperature range above Tc in the UD region of the phase 
diagram but the range ends substantially below the tem- 
perature scale of the onset of pseudogap behavior [7J [8] . 
Alternative models interpret the anomalous properties of 
the pseudogap phase as precursors to the Mott insulator 
at zero doping [T|. Since Mott insulating behavior per se 
is not associated with translational symmetry breaking, 
precursor behavior may not be associated with it either. 

Andreev tunneling has been proposed as a distin- 
guished tool to discriminate between SC pairing fluctu- 
ations and precursor insulating in the pseudogap phase 
in an early paper by Deutscher 9 . (For a review of 
Andreev reflection of pure d-wave superconductor see 
Ref[lO] ) It was pointed out that the voltage (or energy) 
scale in Andreev tunneling experiments on UD cuprates 
in the pseudogap phase [11] was substantially below that 
observed in Giaever or single particle tunneling experi- 



ments. The two voltage scales however were the same in 
the overdoped (OD) region. In a subsequent review of 
the Andreev experiments Deutscher concluded that "the 
balance was tilted somewhat against the preformed pairs 
scenario" [T^]. However, the opposite conclusion, namely 
that higher pseudogap energy scale reflects the pairing 
strength while a second lower scale the SC condensation 
energy, was argued for in a later review by Hufner and 
coworkers which examined many experimental results us- 
ing different techniques 13J. To make progress in this 
debate one needs to move beyond qualitative arguments 
about many individual experiments and on to more ex- 
plicit models which can be used to consistently analyze 
a whole range of experiments. 

In this letter, we perform such solid analysis by using 
the model proposed by Yang, Rice and Zhang (YRZ), 
which has successfully been applied to explain many 
other experiments, to study both Andreev reflection and 
Giaever type tunneling for the UD cuprates. Good 
agreements are achieved, with both the Andreev reflec- 
tion where a small SC gap As is reported, and Giaever 
type tunneling experiment where a large pseudogap is 
reported. Our theory provides a semi-quantitative de- 
scription of the two gap scenario. 

In the YRZ model, a single particle propagator was 
proposed for the pseudogap phase [14]. The YRZ prop- 
agator was inspired by an analysis by Konik, Rice and 
Tsvelik ^^ of ^ 2D array of lightly doped 2-leg Hubbard 
ladders which gave a set of hole pockets with an energy 
gap on flxed lines in fc-space connected by elastic Umk- 
lapp scattering processes. The adaptation of this anal- 
ysis to a 2D Hubbard model for a square lattice lightly 
doped away from half-fllling was influenced by the func- 
tional renormalization group results at weak to moder- 
ate interaction strength of Honerkamp, Salmhofer and 
collaborators [16] and by Zhang et al [17] early anal- 
ysis of Anderson's resonant valence bond (RVB) pro- 
posal (121 HH] . The strong coupling t- J model was ana- 
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FIG. 1: Panel (a): Fermi surface in YRZ model: Blue pocket 
for UD {x = 0.10) and red curve for OD [x = 0.25). The 
width on the pocket (blue) represents spectral weight of quasi- 
particle. Plane (b): DOS in the model for UD and for OD, 
with the gaps denoted by arrows. 




FIG. 2: Schematic plot for a metal-superconductor junction. 
Left: normal metal; Mid, scattering center, and Right: high 
Tc Cu-oxide superconductor. Tn and Ta represent single 
particle and Andreev reflection tunnelings, respectively. 



lyzed using a Gutzwiller renormalized mean field theory 
(RMFT). In this phenomenological approach the single 
particle gap at the antinodes is controlled by the RVB 
gap which truncates the tight-binding Fermi surface into 
4 pockets centered on the nodal directions (see Figflk]). 
The d-wave SC energy gap opens up primarily on these 
Fermi pockets. 

This phenomenological propagator gives a consistent 
description [M] of angle resolved photoemission spec- 
troscopy (ARPES) experiments which followed the evolu- 
tion of the Fermi surface from 4-disconnected Fermi arcs 
centered on the nodal directions in UD to the full Fermi 
surface in OD cuprates [20 . The model has been used 
to explain a range of other spectroscopic measurements 
[21], e.g. ARPES results showing increasing particle- 
hole asymmetry as one moves away from the nodal direc- 
tions [22j, angle integrated photoemission electron spec- 
troscopy (AIPES) experiments measuring the doping, x , 
dependence of the density of states (DOS) [23^ and scan- 
ning tunneling microscopy (STM) measurements of the 
coherent Bogoliubov quasiparticle dispersion at low tem- 
peratures [24 . The YRZ form was also used by Car- 
botte, Nicol and coworkers to successfully describe the 
T and x dependences of a wide range of properties in 
the pseudogap phase including specific heat |25j . opti- 
cal conductivity [26 . London penetration depth |27j and 
symmetry dependent Raman scattering spectra [55| . The 
latter was also analyzed in a similar way by Valenzuela 
and Bascones (29j . 

In the YRZ model [13] , the incoherent part of the single 
particle Green's function contributes a smooth spectral 
background with a tiny real component at low energies, 
and the coherent part reads. 



G°(fc,c^)=.9t/[L.-efc-E°(fc,c^)] 



(1) 



where gt — 2a;/(l + a;) is a renormalization factor [T7lll9j. 
We use the result from the RMFT for the RVB state as 
the "bare" dispersion e^ [SU] . S" is the self energy, which 
is zero for the OD cuprates {x > Xc = 0.2). For the UD 



cuprates, x < Xc 
e^ = — 2t(cosfc3; ^ 
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a;/a;c)(cosfcj; — cosfcj,), with Ao = 0.3. All energies are 
in unit of tg = 0.3eV. Note that e^ = at the antifer- 
romagnetic reduced Brillouin zone boundary, where the 
Umklapp scattering is strongest [16 . Eq[l] predicts four 
Fermi pockets in the pseudogap phase, consistent with 
the recent laser ARPES data [3T]. Fig [Ik] shows one of 
these Fermi pockets. 

Here we consider a d-wave superconducting gap func- 
tion, A^'^ — A^'^ {cos kx — cos ky) for the states around the 
Fermi surface within a small energy shell. Several previ- 
ous work proposed similar SC pairing form of YRZ model 
[H [211 [21]. We choose Ag'^ = O.OSto for UD {x = 0.1), 
and Ag^ = 0.04io for OD {x = 0.25). These gap param- 
eters lead to the SC gap As (maximum | A^'^ | on Fermi 
surface) comparable to those observed in the Andreev 
reflection (~15 and 21 meV, respectively) [11], and to 
those reported in the recent STM data [211 . The Green's 
functions for the SC state in the UD and OD cuprates 
take the following form in Nambu spinor representation. 
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where r. 



/y/'- 



are the Pauli matrices, and the label i = 1, 2 



denotes the lower and higher energy quasiparticle bands 
given by Eq[l| for the UD region with spectral weight 
gtZk,i and dispersion E^^i- Gqj^ has a BCS form. In the 
UD region, the quasiparticle energy of the band i = 2 
is well above the chemical potential, which supports our 
choice A|';2 = 0. The typical profile of DOS of the SC 
state is shown in Fig[lb]. In the UD case, the peaks at 
±90 meV are related to the RVB gap around the antin- 
odes. 

We use the Keldysh formalism to study Andreev re- 
flection of a NS junction, consisting of a normal metal 
lead at the left and a superconductor lead at the right as 



shown in Fig[2]. We approximate the connection area or 
Une interface as a scattering center. The tunneUng of an 
electron from the left(L) to the right (R) via the scatter- 
ing center can be described by the hybridization Hamil- 

tonian, H' = Z]fc.Q=L/fl(^fc4.a'^ + ^■^■) '^ith Ck,a and d 
the electron annihilation operators at lead a = L, R and 
on scattering center, respectively. In our following calcu- 
lations, we assume i^ = i" for simplicity. In an Andreev 
reflection process (for a review see Ref[32]), an electron 
in the normal metal moving along the y-direction to the 
scattering center may be reflected back as a hole due to 
the proximal superconductivity in the scattering center. 
In the tunneling process, kx is conserved. As studied by 
Meir et al. [33] for NN junction, and by Sun et al. [34] 
and Wang et al. [3S] for NS junction, the total current 
can be given by (setting /i = 1), 
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where /" is the Fermi function at lead a. T^^, Ta are 
the diagonal components of the single particle tunneling 
matrix Tjy and the Andreev reflection coefficient, respec- 
tively, s = -l-(— ) corresponds to the electron (hole) chan- 
nel. With the frequency-dependence implicitly implied 
we have, 

TNikx) = G^(fc,)r^(fc,)G^(fc,)r^(fc,) (4) 

TA{kx) - [G^(fc,)]l2[r^(fc.)]22[G^(fc.)]2l[r^(fc.)]ll 

where Gc refers to the retarded or advanced Green's 
function, Gc{kx) = l/[w - ^^{kx) - ^'^(kx)] is the 
Green's function in the scattering center [36]. The self 
energies ^"{kx) — {t°')'^TzG"{kx)Tz and lifetime broad- 
ening V'ikx) = -i[S"'"(A;,) - S«''-(fc,)], with G"(fc,) 
the Green's functions on the metal side or on the su- 
perconductor side given by Eq [l|2j at their edges to the 
scattering center. Spin rotational invariance leads to 
T^'^{uj) = T^~{-uj) and TA(a;) = TAi-uj). 

The conductance at zero temperature is given by 



asieV) 
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[U'{kx, seV) + 2TA{kx, seV)] .(5) 



In our calculations, we consider the tunneling along 
an antinodal direction of the Cu02 plane, and as- 
sume a parabolic dispersion on the normal metal, e^ ~ 
{k^ — kp)/2m, setting kp — ii 12, m = n 12 in unit of 
l/io. 

We now discuss the tunneling conductances in our 
model and compare them with experiment. There are 
two parameter regions in the tunneling. At large i^/i, 
we have transparent tunneling, Andreev reflection is 
dominant within the SC gap A^. At small i^/i, the An- 
dreev reflection is strongly suppressed, and single particle 
tunneling dominates. 
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FIG. 3: Tunneling conductance for OD cuprate (x = 0.25). 
Panel (a) transparent limit. Panel (b) Giaever limit, a^ ipn) 
for NS (NN) junction. O;, is the contribution to Os from 
Andreev reflection. The SC gap (As ~ 21 meV) is indicated 
with green lines. Black dots in Panel (b) is STM dl /dV for 
overdoped Bi2212 (68K) [32]. (broadening <5 = 0.005. t^ = 
g^' t^ — 10, 0.02 for the two limits, respectively.) 



We first discuss the OD case, which is similar to the 
conventional BCS superconductor. In this case, the nor- 
mal state is a metal with a full Fermi surface shown in 
Fig[I^]. The RVB gap vanishes, leaving only a SC gap 
As ■ The calculated conductance as in a typical transpar- 
ent region is shown in Fig[3k]. The conductance in the 
single particle tunneling region is shown in Fig[3L)] in a 
good agreement with STM data [37] , where the contribu- 
tion from the Andreev reflection is essentially vanishing. 

For the UD cuprates, there are two energy scales, a 
larger RVB gap or pseudogap, and a smaller SC gap A^. 
We expect two distinct energy scales in the tunneling 
conductance, as in the transparent region is plotted in 
Figjik]. The contribution from the Andreev reflection is 
only substantial within the SC gap on the remnant re- 
gion of the Fermi surface, as shows a clear peak-edge 
feature at the SC gap energy As. Our result is in good 
agreement with the reported Andreev tunneling experi- 
ment |11| , which is reproduced here for comparison. In a 
single gap scenario, it is not clear that the relative contri- 
butions of Andreev and single particle processes to tun- 
neling, into the ordered SC state should change between 
OD and UD. 

The voltage-dependent conductance in the single parti- 
cle tunneling region is plotted in Fig[4b] for UD samples. 
In this case, the Andreev reflection is substantially sup- 
pressed, so that voltage-dependence of as is similar to 
that of DOS (Fig[l]3]). In both as and DOS in the SC 
phase, there are two energy scales. The lower energy peak 
is associated with the SC gap Ag, and the higher energy 
peak with the RVB gap. The overall profile we obtained 
is very close to the recent STM experimental data on Bi- 
2212 [38]. Note that the tunneling conductance we study 
here does not include a strong electron-hole asymmetry 
in the spectral weight associated with the asymmetry of 
injecting an electron or a hole in the strongly correlated 
systems, as examined by Anderson and Ong ,39) and 
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FIG. 4: Conductance for UD cuprate {x = 0.1). Panel (a) 
transparent limit. Red curve is the NS conductance as, a^ in 
blue denotes the Andreev reflection contribution, black curve 
is experimental Andreev reflection data on UD YBCO sam- 
ples [11]. The SC gap (As ~ 15 meV) is indicated with green 
lines. Panel (b) Giaever tunneling limit. Blue curve is a„ 
for NN junction. Black curve is dl/dV curve for UD Bi-2212 
(Tc = 45 K) observed in STM 138]. The two energy gaps 
are marked by arrows. The same parameters as OD shown in 
Figis] are used. 



these antinodal regions through Cooper channel scatter- 
ing processes, the corresponding Andreev signal is found 
strongly suppressed. Since the tunnehng experiments, 
per se, do not address the origin of this two gap behav- 
ior, from broader viewpoint, our calculations support all 
models in which the pseudogap at underdoping is due to 
the partial truncation of the Fermi surface through an 
insulating gap in the antinodal regions, but not due to 
preformed Cooper pairs. 
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by Randeria et al. [40]. We argue that the relatively 
low energy tunneling spectral is weakly affected by this 
strong asymmetry [3D] so that our model calculations of 
the tunneling conductance remain valid in the relevant 
energy region. 

One point worth noting concerns the asymmetry in 
the single particle tunneling spectra between positive and 
negative voltages. This is much more pronounced in the 
DOS calculations based on the YRZ propagator than in 
the tunneling experiments. The asymmetry in the YRZ 
DOS comes from the Dirac point at positive energies in 
the lower quasiparticle band. Interestingly a linear de- 
pendence of the DOS at the chemical potential on the 
hole density was reported in AIPES ^23^, which is con- 
sistent with the approaching to a Dirac point as the hole 
density is reduced. The quasiparticle DOS calculations 
do not include the effects of lifetime broadening of the co- 
herent quasiparticle spectra. This could be important in 
view of the evidence for strong inelastic scattering pro- 
cesses. These could act to smear out signs of a Dirac 
point at finite energies above the chemical potential and 
so reduce the asymmetry in the DOS. This discrepancy 
requires further study. 

In summary, we have theoretically studied the tunnel- 
ing spectroscopy of the underdoped cuprates. The single 
particle pseudogap in our theory appears as a precursor 
to the Mott insulator at zero doping. Our theory semi- 
quantitatively explains the two gap scenario at under- 
doping revealed in tunneling experiments: a small energy 
gap associated with Andreev refiection in the transpar- 
ent limit and a large gap associated with single particle 
tunneling, which are in contrast with the same scales 
at overdoping [S] [T^] . Although the onset of long range 
superconducting order on the remaining Fermi pockets 
(or arcs) induces a reduced pairing amplitude also in 
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